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This  paper  presents  a  procedure  to  evaluate  the  optimal  element  sizing  of  hybrid  power  systems.  In  order 
to  generalize  the  problem,  this  work  exploits  the  “energy  hub”  formulation  previously  presented  in  the 
literature,  defining  an  energy  hub  as  an  interface  among  energy  producers,  consumers  and  the  transporta¬ 
tion  infrastructure.  The  resulting  optimization  minimizes  an  objective  function  which  is  based  on  costs 
and  efficiencies  of  the  system  elements,  while  taking  into  account  the  hub  model,  energy  and  power 
constraints  and  estimated  operational  conditions,  such  as  energy  prices,  input  power  flow  availability 
and  output  energy  demand.  The  resulting  optimal  architecture  also  constitutes  a  framework  for  further 
real-time  control  designs. 

Moreover,  an  example  of  a  hybrid  storage  system  is  considered.  In  particular,  the  architecture  of  a 
hybrid  plant  incorporating  a  wind  generator,  batteries  and  intermediate  hydrogen  storage  is  optimized, 
based  on  real  wind  data  and  averaged  residential  demands,  also  taking  into  account  possible  estimation 
errors.  The  hydrogen  system  integrates  an  electrolyzer,  a  fuel  cell  stack  and  hydrogen  tanks.  The  resulting 
optimal  cost  of  such  hybrid  power  plant  is  compared  with  the  equivalent  hydrogen-only  and  battery-only 
systems,  showing  improvements  in  investment  costs  of  almost  30%  in  the  worst  case. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  energy  infrastructures  of  today  are  about  to  undergo  a  pro¬ 
found  change:  fossil  fuel  prices  are  raising  every  year  while,  at  the 
same  time,  energy  demand  increases  in  every  country.  Moreover, 
the  aim  to  reduce  greenhouse  gas  emissions  is  moving  its  attention 
to  more  environmentally-friendly  and  sustainable  energy  sources. 
With  an  increased  utilization  of  small  distributed  energy  resources 
for  generation  of  electricity  and  heat  [1  ],  renewable  energy  genera¬ 
tion  will  constitute  an  important  part  of  the  overall  energy  scenario 
in  the  coming  years. 

One  of  the  main  problems  associated  with  these  kind  of  sys¬ 
tems  is  the  reliability  and  quality  of  the  power  supply.  As  a  matter 
of  fact,  since  the  renewable  source  is  intermittent,  unpredictable 
fluctuations  may  appear  in  power  output  [2].  Also,  electrical  gen¬ 
eration  from  renewable  sources  is  not  subject  to  demand,  which 
creates  imbalance  in  the  system.  One  way  to  overcome  this  problem 
is  by  including  intermediate  storage,  such  as  batteries,  water  pump¬ 
ing,  super-capacitors,  compressed  air,  fly  wheels,  superconducting 
magnetic  energy  storages,  etc.  [3].  Among  the  most  promising 
storage  technologies  are  those  based  on  hydrogen  production  and 
utilization,  which  is  expected  to  be  used  for  very  different  appli¬ 
cations  [4,5]  as  they  constitute  some  interesting  advantages  in 
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terms  of  cost,  autonomy,  power  range  and  environmental  effects 

[6]. 

However,  hybrid  energy  storage  systems  increase  the  com¬ 
plexity  of  the  overall  power  plant,  the  control  design  having  an 
important  effect  on  system  performance.  Thus,  there  are  a  num¬ 
ber  of  controllers  available  in  the  literature,  such  as  those  based 
on  heuristic  rules  and  trial-and-error  techniques  [7,8].  Fuzzy  logic 
approaches  [9,10]  are  equivalent  to  those  based  on  heuristic  rules 
in  the  sense  that  they  rely  on  system  knowledge  to  obtain  the 
‘best’  intuitive  power  management.  Nonetheless,  other  approaches 
based  on  on-line  optimization  can  be  found,  resulting  in  a  more 
re-usable  and  rigorous  design  process,  so  that  the  final  algorithm 
achieves  a  guaranteed  optimum  level  [11].  Along  these  lines,  an 
on-line  optimization  to  minimize  the  hydrogen  consumption  for 
residential  hybrid  power  plants  was  presented  in  [12,13].  As  for 
renewable  sources,  the  intermittency  of  the  available  power  also 
has  a  great  impact  on  system  performance.  Although  not  being 
suitable  for  real-time  control  as  the  designs  cited,  there  are  some 
control  algorithms  based  on  prior  knowledge  of  future  conditions 
(such  as  wind  speed  data)  which  are  useful  as  a  basis  of  comparison 
for  the  evaluation  of  real-time  control  strategy  quality  [14]. 

As  well  as  the  control  design,  it  is  very  important  that  compo¬ 
nent  sizing  be  taken  into  account  in  order  to  reduce  installation 
investment  costs  and  to  achieve  good  overall  performance,  which 
is  the  main  objective  of  this  work.  However,  very  few  papers  have 
addressed  this  issue.  To  this  end,  [15]  discusses  the  best  coupling 
methods  for  conventional  storage  batteries  with  hydrogen  energy 
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Fig.  2.  Energy  hub  basic  elements:  converter  (left)  and  storage  (right). 


Fig.  1.  General  energy  hub  diagram. 


storage  which  includes  an  electrolyzer,  hydrogen  storage  tank,  and 
a  fuel  cell.  The  resulting  study  shows  that  if  multiple  energy  storage 
devices  with  complementary  performance  characteristics  are  used 
together,  the  resulting  hybrid  system  can  dramatically  reduce  the 
cost  of  energy  storage  over  single  storage  systems.  Although  [15] 
is  a  good  study  about  element  sizing,  it  is  only  aimed  for  a  certain 
power  system,  which  makes  very  difficult  to  analyze  other  differ¬ 
ent  system  layouts.  In  contrast,  this  paper  utilizes  a  very  general 
formulation,  which  can  be  used  to  model  and  study  any  possible 
layout. 

In  particular,  the  general  sizing  problem  presented  in  this  paper 
is  a  revision  of  the  work  presented  in  [16]  and  other  related  papers 
by  the  same  authors  such  as  [14,17-19],  which  proposes  a  general 
mathematical  formulation  to  model  hybrid  power  plants,  introduc¬ 
ing  the  so-called  “energy  hubs”.  An  energy  hub  is  defined  as  an 
interface  among  energy  producers,  consumers  and  the  transporta¬ 
tion  infrastructure.  Specifically,  [16]  finds  the  optimal  hub  layout 
of  a  power  system  but  does  not  discuss  the  sizing  problem,  while 
the  formulation  changes  introduced  herein  are  aimed  to  determine 
the  optimal  hub  size  of  a  certain  layout,  thus  completing  the  men¬ 
tioned  previous  work.  Also,  as  other  of  this  work  contributions,  a 
particular  hybrid  hydrogen-based  power  system  is  optimally  sized 
based  on  real  wind  and  consumption  data.  Such  system  is  derived 
from  the  general  sizing  optimization  framework  proposed. 

In  the  following  section,  energy  hub  concept  and  mathematical 
formulation  are  briefly  outlined,  as  there  is  an  extensive  literature 
by  the  corresponding  authors  describing  them.  Section  3  proposes 
a  general  cost  function  to  minimize  component  sizing  based  on 
costs  and  efficiencies.  Section  4  applies  the  general  optimization 
scenario  to  a  wind  generator/hydrogen/batteries  power  plant,  also 
discussing  the  results  obtained.  Lastly  Section  5  is  dedicated  to  the 
concluding  remarks. 

2.  Energy  hub  concept  and  formulation 

As  an  increased  utilization  of  distributed  generation  technolo¬ 
gies  will  characterize  future  energy  systems,  terms  like  “multiple 
energy  carrier  systems”  [20]  and  “hybrid  energy  systems”  [21  ]  have 
become  the  norm  when  referring  to  systems  including  various 
forms  of  energy.  In  this  way,  as  noted  in  [19],  there  are  a  number 
of  approaches  to  formulate  these  kind  of  systems,  such  as  “energy- 
services  supply  systems”  [22],  “basic  units”  [23],  “microgrids”  [24] 
and  the  so-called  “hybrid  energy  hubs”  [25]. 

The  latter  formulation  is  adopted  herein,  which  is  extensively 
described  in  the  PhD  thesis  [19]  and  related  publications.  Accord¬ 
ing  to  this  formulation,  energy  hubs  are  defined  as  interfaces  among 
energy  producers,  consumers,  and  the  transportation  infrastruc¬ 
ture  (see  Fig.  1,  where  P;  are  power  inputs  and  Lj  power  outputs), 
and  contain  three  basic  elements:  direct  connections,  converters 
and  storage  (see  Fig.  2,  with  Q/<  being  the  power  exchange,  Q/^  the 
internal  power  and  Ek  the  stored  energy). 

Converters  link  inputs  and  outputs  through  coupling  factors  qj, 
which  can  be  considered  to  be  the  converter’s  steady-state  energy 


efficiency,  expressed  as: 


Lj  —  ^ijPi 


(1) 


Considering  all  the  energy  hub  inputs  P  and  outputs  L,  the  following 
converter  coupling  matrix  C  results: 


(2) 


-ir 

‘Cl,l 

cf,r 
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-U- 

-cu  ••• 

Ci,j  - 
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P 

As  the  input  flow  P*  can  be  distributed  among  various  converter 
devices  (see  Fig.  3),  dispatch  factors  vi  n  specify  how  much  of  the 
input  power  Pt  flows  into  the  converter  n: 

Pi,n  —  nPj  (3) 

Conservation  of  power  also  introduces  the  constraints 
0  <  vUn  <  1  Vi,  Vn  (4a) 

)Tv„n  =  1  Vi  (4b) 


With  respect  to  storage,  power  exchange  Q/<  and  stored  energy  Ek 
are  linked  through  the  equation: 


(5) 


ek  being  the  efficiency  of  the  charge/discharge  storage  interfaces, 
expressed  as 


ek 


{> 


=  i  et 


tel 


if Q/<  >  0  (charging/ standby) 
else  (discharging) 


(6) 


When  storage  elements  exist,  power  conservation  leads  to  the  fol¬ 
lowing,  depending  on  which  side  of  the  converter  the  storage  is 
located  (see  Fig.  4): 


Pi  =  Pi-Qi 


Lj  ~  Lj  +  LAj 


Adding  the  storage  to  the  hub  Eq.  (2)  leads  to: 
[L  +  M]  =  C[P  —  Q] 


(7a) 

(7b) 

(8) 
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Fig.  4.  Converter  with  storage  at  the  input  and  the  output  sides. 


Assuming  a  constant  converter  coupling  matrix  C  and  applying 
superposition,  the  equivalent  storage  flows  are: 

Meq  =  CQ+  M  (9) 

Rewriting  (8)  in  a  more  condensed  form, 

L  =  C  P  -  Meq  (10) 

Defining  the  storage  coupling  matrix  S  to  describe  how  changes 
of  the  storage  energy  derivatives  affect  the  hub  output  flows,  the 
equivalent  storage  power  flows  Meq  can  be  stated  as 


-JW^q- 

"$1,1 

•••  S1  ,k~ 

_S1  ,k 

$k,k  _ 

Meq 

S 

Summarizing  all  the  previous  equations,  the  complete  hub  energy 
model  would  be: 

L  =  CP  -  SE  (12) 

3.  Optimal  hub  size 

Optimal  hub  design  can  be  divided  into  two  different  steps:  opti¬ 
mal  hub  sizing  and  hub  control  design  (see  Fig.  5 ).  Most  of  the  papers 
in  the  literature,  as  mentioned  in  Section  1,  are  dedicated  solely  to 
the  controller  design,  not  addressing  the  sizing  problem.  The  better 
the  controller  design,  the  better  the  performance  of  a  given  system. 
However,  hub  sizing  and  control  design  are  not  independent  from 
one  another.  As  a  matter  of  fact,  the  performance  of  the  overall  sys¬ 
tem  not  only  depends  on  the  quality  of  the  controller  but  also  on 
the  hub  characteristics. 

Optimal  hub  sizing  for  any  given  hub  layout  entails  optimiza¬ 
tion  of  converter  and  storage  element  sizes.  To  that  end,  cost  and 


OPTIMAL  ENERGY  HUB  DESIGN 

Fig.  5.  Optimal  energy  hub  design  steps. 


efficiencies  associated  with  each  component,  as  well  as  the  esti¬ 
mated  working  conditions  of  the  hub  (such  as  energy  prices,  input 
energy  flows  availability,  output  power  demand,  etc.)  have  to  be 
taken  into  account.  Such  estimated  working  conditions  are  usually 
based  on  historical  data.  However,  since  the  optimization  is  done 
on  estimated  power  inputs  and  outputs,  estimation  errors  are  likely 
to  occur.  Thus,  also  the  sensitivity  of  the  computed  optimal  size 
solution  to  such  errors  should  be  taken  into  account  considering 
a  worst-case  power  availability  (input)  and  consumption  (output), 
which  can  be  done  computing  the  optimization  for  different  esti¬ 
mated  conditions  data  sets  and  selecting  the  most  conservative  hub 
sizing. 

Given  the  optimal  hub  sizing,  a  suitable  optimization  problem 
minimizing  a  determined  objective  function  and  supposing  a  set 
of  system  operational  conditions  would  then  represent  the  basis  of 
comparison  for  the  evaluation  of  real-time  control  strategy  quality. 
This  type  of  optimization  problem  is  referred  to  as  “optimal  power 
dispatch”  [14](also  shown  in  Fig.  5). 

The  problem  presented  by  optimal  hub  sizing,  which  is  the 
objective  of  this  work,  can  be  basically  expressed  with  three  rela¬ 
tions:  physical  laws  representing  the  hub,  its  technical  limitations, 
and  an  objective  function  which  accounts  for  the  minimization  of 
the  system  investment  cost.  Then,  the  optimization  is  stated  as  a 
multi-period  nonlinear  constrained  problem  including  an  objective 
function,  equality  and  inequality  constraints. 

The  physical  laws  representing  the  energy  hub  are  modeled  by 
the  equality  constraints  presented  in  Section  2.  Extending  that  for¬ 
mulation  to  consider  multiple  time  periods,  the  hub  model  equation 
would  be: 

L(t)  =  c(0  p(t)  _  $(0  £(f)  vf  (13) 

Inequality  constraints  correspond  to  the  technical  limitations  of 
the  converter  and  storage  elements.  Eq.  (14a)  expresses  power 
conversion  limitations  of  the  converters.  Eqs.  (14b)  and  (14c)  cor¬ 
respond  to  charging  and  discharging  power  limits  of  the  storage 
interfaces,  while  (14d)  considers  the  energy  capacity  limits  of  the 
storage  devices.  The  last  inequality  (14e)  is  also  included  so  that 
stored  energy  at  the  end  of  the  optimization  period  Nt  is  equal  to 
or  greater  than  the  initial  amount,  in  order  to  ensure  sustainable 
storage  utilization,  where  Nt  is  the  last  considered  period  of  the 
optimization. 


<PUn  Vt,  Vi,  Vn 

(14a) 

Q,  <  Q,- 0  <  Qj  Vt,  Vi 

(14b) 

Mj  <  m]°  <  Mj  Vt,  Vj 

(14c) 

£fc  <  Ef  <  Ek  Vt,  Vk 

(14d) 

E(°)  <  E(Nt)  yk 

(14e) 

The  objective  function  T  depends  on  converter  and  storage  ele¬ 
ment  sizes  Pi  n,  Ek,  Qj,  Q..,  Mj  and  Mj,  which  are  related  to  the 
maximum  power  conversion  and  storage  capacities.  Notice  that  Q.. 
and  Mj  are  the  maximum  discharging  rate  of  the  storage  interfaces, 

as  Q jc\  Mf  <  0  during  discharges  and  thus  Q..,  Mj  <  0.  Converter 
and  storage  elements  sizes  appear  in  the  cost  function  multiplied 
by  coefficients  Cp.,  c c^.,  Cq.,  and  Cm.,  which  account  for  the 
investment  cost  of  hub  elements  depending  on  their  size.  Thus,  the 
solution  of  the  optimization  problem  provides  the  optimal  sizes 
of  the  hub  elements,  minimizing  the  total  investment  cost  while 
taking  into  account  the  hub  model  equations  and  its  technical  lim¬ 
itations. 
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Considering  a  quadratic  function,  the  objective  can  be  expressed 
as: 

?  =  + DcqA2 + 

i  k  i 

+  J2(CMjMf+CMM ?)  (15) 

j 

The  hub  size  optimization  problem  can  finally  be  stated  as: 


Minimize  objective  function 
subject  to  energy  hub  model 

energy  and  power  constraints 

When  the  objective  function  is  convex  and  the  constraints  are 
expressed  as  linear  equations,  the  global  optimum  can  be  found 
utilizing  numerical  methods,  as  the  solution  space  is  convex. 

This  approach  can  be  considered  as  one  of  the  main  contribu¬ 
tions  of  this  work,  as  it  presents  a  general  optimization  framework 
for  element  sizing,  completing  the  work  initiated  in  [16]  and 
improving  [15]  as  it  only  studied  a  certain  system,  not  a  general 
framework  for  any  system  layout  as  here.  Specifically,  the  formula¬ 
tion  changes  with  respect  to  [16]  are  related  to  the  element  sizes, 
which  are  considered  here  as  optimization  variables  appearing  in 
the  objective  function  and  also  in  the  constraint  equation  set,  as  an 
innovative  approach  to  the  sizing  problem. 


Fig.  7.  Corresponding  energy  hub  of  a  hybrid  energy  storage  system. 


Table  1 

Hybrid  energy  storage  system  efficiencies. 

Hub  element 

Efficiency 

Electrolyzer 

tie  =  0-74 

Fuel  cell 

Ore  =  0.47 

Battery  charging 

Ob,  ch  =  °-7 

Battery  discharging 

0b,  dis  =  0-9 

EW  =  [£W 


Efl 


(18) 


Also,  following  the  aforementioned  notation,  converter  coupling 
matrix  C(t)  and  storage  coupling  matrix  S(t)  are  stated  as: 


4.  Application 


C(t)  =  [l] 


(19) 


Considering  the  system  shown  in  Fig.  6,  the  primary  energy 
source  is  a  wind  generator,  which  is  connected  to  a  residential  load 
(Lr).  The  electricity  produced  via  wind  (w)  can  be  delivered  to  the 
load  and/or  be  diverted  to  an  electrolyzer  (E)  and  batteries  (B).  The 
energy  consumed  by  the  electrolyzer  (Qf )  is  used  to  produce  hydro¬ 
gen,  which  is  stored  in  the  tanks  placed  in  the  hydrogen  line  (Eh2  ). 
The  fuel  cell  stack  (FC),  fed  by  those  tanks,  can  produce  electric¬ 
ity  (Qfc)-  Similarly,  the  batteries  can  be  charged  (Qb, chX  storing  the 
energy  (EB),  and  discharged  {Qb, dis\  thus  complementing  the  total 
power  supplied  to  the  load. 

Deriving  this  specific  case  from  the  general  problem,  and  assum¬ 
ing  a  certain  set  of  operational  conditions,  the  optimal  hub  sizing 
for  the  proposed  system  is  calculated.  To  that  end,  the  optimization 
problem  is  formulated  as  in  the  previous  Section  3. 

4.1.  Energy  hub  model 

Model  equations  are  based  on  the  notation  presented  in  Section 
2.  This  way,  the  specific  energy  hub  is  illustrated  in  Fig.  7.  Input,  out¬ 
put  and  storage  energy  derivative  vectors  for  multiple  time  periods, 
can  be  defined  as 


II 

(16) 

II 

(17) 

Eh, 

Fig.  6.  Hybrid  energy  storage  system. 


S(t)  =  [l/e£  \/ef\  (20) 

where  and  ef  are  the  storage  interface  efficiencies,  the 
electrolyzer  being  the  ‘charging’  interface  and  the  fuel  cell  the  ‘dis¬ 
charging’  interface  for  the  hydrogen  line.  Also,  notice  that  different 
battery  charging  and  discharging  efficiencies  are  considered  (see 
Table  1),  resulting  in  the  following  relations: 


,(0  _  f  >1E 

if  ojfj  >  0 

(electrolyzer) 

" 2  {  1/*|FC 

else 

(fuel  cell) 

,(0  _  /  'iB.ch 

ifQ^0  >  0 

(battery  charging) 

B  \  V'lB.dis 

else 

(battery  discharging) 

(21a) 

(21b) 


with  the  power  exchanges  Q^,  Q^f)  and  storage  energy  derivatives 
)  expressed  as: 


^=eH2<)-4"1)Qlt_1)  (22a) 

4t)=eB)(^t)-eB_1)(^t_1)  (22b) 


4.2.  Energy  and  power  constraints 

Technical  limitations  are  modeled  as  in  (14).  In  this  way,  input 
power  limits,  storage  interfaces  power  exchange  capacities  and 
stored  energy  limitations  are  evaluated  next. 

With  respect  to  the  input,  the  constraint  vector  stated  in  (23a) 
incorporates  the  wind  generator  model,  taking  into  account  that 
if  the  wind  becomes  too  strong,  the  wind  generator  needs  to  be 
switched  off  and  no  wind  in-feed  is  apparent  to  be  injected  into 
the  system  (see  Fig.  8).  Introducing  the  function  cw,  Ws(t)  being  the 
actual  wind  speed  and  being  the  power  generated  by  the  wind 
generator,  the  input  power  constraint  results: 

[0]  <  [P(f]  <  [cwWs(t)] 


(23a) 
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Fig.  8.  Normalized  wind  generator  curve. 


Cw 


0 

l/10-l/(2Wf) 

1/Wf 

0 


if  wjf)  <  5 
if  5  <  Ws(t)  <  15 
if  15  <  wf  <  25 
if  wf  >  25 


(23b) 


Power  storage  exchange  is  also  limited  by  the  maximum  power  that 
can  be  provided  by  the  storage  interfaces: 


-Qfc 

< 

ro^i 

< 

Qe 
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Lq^J 

Qb, ch 

(24) 


whereas  for  the  hydrogen  line,  Qf  and  Qfc  represent  the  maximum 
capacities  of  electrolyzer  and  fuel  cell,  respectively.  Concerning  the 
batteries,  QB  Ch  and  Q^^is  are  the  limit  charging/discharging  rates. 
Notice  that  these  rates  are  usually  a  function  of  total  battery  size  £, 
assuming  here  that  Qgch  =  0.2 E  and  Qb, dis  =  IE- 

Maximum  stored  energy  depends  on  the  size  of  the  hydrogen 
tanks  Eh2  and  the  batteries  EB.  Due  to  technical  constraints,  the  bat¬ 
teries  should  never  be  totally  drained  nor  fully  charged;  they  should 
always  be  in  a  partially  charged  state.  Taking  these  considerations 
into  account  and  forcing  a  safe  charge  level  of  (0.2 EB,  0.9 EB),  the 
constraint  can  be  expressed  as: 
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Finally,  the  following  constraint  to  verify  sustainable  energy  storage 
is  also  introduced: 


< 


(26) 


4.3.  Objective  function 

Moving  from  the  general  (15)  to  the  specific,  the  objective  func¬ 
tion  would  be: 

Jr=  c~pJI  +  %2  Efi2  +  %Eb  +  cQfCQfc  +  cQe^e  (27) 

where  the  cost  coefficients  are  based  on  the  investment  costs  pre¬ 
sented  in  [15](Table  2). 


Table  2 

Hybrid  energy  storage  element  costs. 


4.4.  Operational  conditions 

As  the  optimal  hub  architecture  design  is  based  on  estimated 
operational  conditions,  the  more  precise  the  utilized  data  are,  the 
more  accurate  are  the  architectural  results.  As  for  the  power  input, 
three  different  wind  power  normalized  data  sets  were  considered 
(see  Fig.  9),  which  correspond  to  data  recorded  during  different 
periods.  Concerning  the  load  Lr,  the  data  used  is  shown  in  Fig.  10, 


Set  (a)  Available  Wind  Power  (Normalized) 


Fig.  10.  Residential  sector  average  daily  loads. 
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Table  3 

Size  optimization  results  for  the  hybrid  storage  system. 


which  represents  the  average  daily  load  for  the  residential  sector  in 
Spain  [26].  The  sampling  time  for  all  the  data  sets  is  1  h. 

4.5.  Optimization  results 

As  the  objective  function  is  convex  and  the  constraints  are 
expressed  as  linear  equations,  the  global  optimum  can  be  found  uti¬ 
lizing  numerical  methods,  as  the  solution  space  is  convex.  Thus,  the 
optimization  was  implemented  in  MATLAB  using  the  commercial 
solver  “CPLEX”,  resulting  in  a  Mixed  Integer  Quadratic  Program¬ 
ming  (MIQP).  Also,  the  optimization  was  performed  for  the  three 
power  input  data  sets  shown  in  Fig.  9  and  the  load  demand  shown 


in  Fig.  10,  in  order  to  size  the  system  in  the  worst  possible  wind 
availability  scenario. 

The  results  obtained  confirm  that  the  third  scenario  (data  set 
(c)  in  Fig.  9)  is  the  most  disadvantageous,  mainly  due  to  the  fact 
that  there  is  less  wind  power  available.  Specifically,  the  investment 
cost  of  such  scenario  is  5.9%  higher  than  the  investment  cost  for 
the  most  advantageous  conditions,  which  correspond  to  data  set 
(a).  The  optimal  size  of  the  hub  elements  and  the  investment  cost 
of  the  installation  for  the  three  data  sets  considered  are  shown  in 
Table  3. 

The  hybrid  storage  system  sized  for  the  most  disadvantageous 
wind  conditions  was  simulated  for  a  different  wind  speed  data  set, 
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Fig.  11.  Optimal  power  and  energy  management  for  a  hybrid  storage  power  system. 


A.J.  del  Real  et  al.  /  Journal  of  Power  Sources  193  (2009)  315-321 


321 


Table  4 

Sensibility  analysis. 

Estimation  error 

Oversizing  requirement 

Peak  load  demand  +1 W 

Battery  +38.4Wh 

No-wind  during  +1  h 

Electrolyzer  +6.47  W 

H2  tank +382.97  Wh 

Table  5 

Hybrid  system  and  other  layouts  comparison. 


Hydrogen-only 

Battery-only 

Equipment 

Cost  (size) 

Cost  (size) 

Wind  generator 

$8600 (4300 W) 

$12500 (6250 W) 

Electrolyzer 

$665  (350  W) 

- 

Fuel  cell 

$1375  (550  W) 

- 

H2  tank 

$3681  (122.695  kWh) 

- 

Batteries 

- 

$3651  (18.255  kWh) 

Total  cost 

$14321 

$16151 

Cost  increment 

11.91% 

28.02% 

as  can  be  seen  in  Fig.  11 ,  which  also  shows  batteries,  electrolyzer  and 
fuel  cell  utilization.  In  particular,  the  fuel  cell  is  used  as  a  base  power 
supplier,  while  batteries  are  utilized  to  deliver  the  power  peaks. 
As  expected,  the  hybrid  storage  system  studied  combines  the  best 
characteristics  of  both  energy  storage  devices.  It  is  also  remarkable 
that  the  system  can  supply  all  the  power  residential  demand  when 
simulated  with  the  wind  speed  data  shown  in  Fig.  11,  as  such  wind 
speeds  resulted  to  be  more  advantageous  than  the  wind  conditions 
used  to  size  the  system  (set  data  (c)  in  Fig.  9).  The  case  that  the  wind 
conditions  were  worse  than  the  conditions  which  the  system  was 
designed  for  is  discussed  in  next  Section  4.6. 

4.6.  Sensibility  analysis 

Although  being  done  for  the  worst  possible  scenario,  estima¬ 
tion  errors  are  likely  to  occur,  i.e.,  if  the  real  conditions  are  even 
worse  than  the  conditions  considered  herein.  Thus,  a  sensibility 
analysis  was  performed  to  take  into  account  these  errors  and  its 
associated  costs.  In  particular,  the  analysis  was  realized  consider¬ 
ing  two  error  estimations:  increments  in  the  load  demand  during 
peaks  and  longer  periods  without  wind  available.  During  such  sit¬ 
uations,  the  optimally  sized  system  would  not  be  able  to  satisfy  the 
power  demands.  One  possibility  to  overcome  these  situations  is  to 
oversize  the  system. 

Specifically,  as  the  simulations  (see  Fig.  11)  show  that  batteries 
are  utilized  to  deliver  the  power  peaks,  the  battery  size  should  be 
incremented  in  order  to  take  into  account  possible  peak  load  esti¬ 
mation  errors.  In  order  to  consider  wind  speed  prediction  errors,  the 
optimization  tests  showed  that  the  electrolyzer  and  the  hydrogen 
tank  should  be  oversized  in  order  to  accumulate  more  hydrogen, 
which  would  be  available  during  periods  without  power  from  wind. 
The  detailed  analysis  results  are  shown  in  Table  4,  taking  into 
account  peak  load  estimation  errors  of  1  W  and  no-wind  extra 
periods  of  1  h. 

4.7.  Different  layouts  comparison 

Finally,  the  optimization  was  done  for  two  other  possible  lay¬ 
outs:  hydrogen-only  and  battery-only  storage,  based  on  the  wind 
power  availability  data  set  (c)  shown  in  Fig.  9  and  the  load  demand 
presented  in  Fig.  10.  In  particular,  hydrogen-only  storage  cost  is 
11.91%  higher  than  the  hybrid  plant,  the  battery-only  choice  being 


28.02%  more  expensive  than  such  hybrid  system  (see  Table  5  for 
detailed  information).  In  a  hydrogen-only  choice,  the  fuel  cell  size 
has  to  be  increased  in  order  to  create  the  power  peaks,  which  results 
in  a  cost  increase  due  to  the  high  cost  of  the  equipment.  On  the  other 
hand,  battery-only  storage  requires  a  large  total  energy  capacity, 
which  is  costly  too.  Also,  as  the  hybrid  system  can  manage  power  in 
a  flexible  and  efficient  way,  wind  generator  for  such  hybrid  system 
can  be  smaller  than  the  wind  generator  needed  by  the  two  other 
layouts. 

5.  Concluding  remarks 

In  this  paper,  an  optimization  strategy  for  sizing  hybrid  power 
systems  is  presented,  improving  the  mathematical  formulation 
based  on  the  “energy  hub”  concept  described  in  previous  litera¬ 
ture.  The  optimization  procedure  was  applied  to  a  hybrid  power 
plant  incorporating  a  wind  generator,  conventional  batteries  and  a 
hydrogen  storage  system  comprised  of  a  fuel  cell,  an  electrolyzer 
and  hydrogen  tanks.  Also,  a  sensibility  analysis  was  performed  to 
take  into  account  possible  estimation  errors.  Finally,  the  optimal 
architecture  was  compared  to  other  possible  layouts,  resulting  in 
almost  a  30%  improvement  among  the  possible  system  layouts  in 
terms  of  cost  effectiveness,  showing  the  flexibility  and  usefulness 
of  the  improved  formulation  proposed  herein. 
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